Thermodynamics of second phase conductive filaments 
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We present a theory of second phase conductive filaments in phase transformable systems; appli- 
cations include threshold switches, phase change memory, and shunting in thin film structures. We 
show that the average filament parameters can be described thermodynamically. In agreement with 
the published data, the predicted filament current voltage characteristics exhibit negative differential 
resistance vanishing at high currents where the current density becomes a bulk material property. 
Our description is extendible to filament transients and allows for efficient numerical simulation. 
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A general observation common to multiple modern 
technologies is that thin-film structures can drastically 
decrease their transversal electric resistance by forming 
conductive filaments (CF) in response to electric bias. 
CF can be either stable, as in phase change memory 
and dielectric oxides after hard breakdown, or unstable 
(disappearing after the bias is removed), as in threshold 
switches [l[ . From a practical perspective, CF can cause 
detrimental shunting and loss of functionality in devices 
such as thin-film photovoltaics and thin oxides of elec- 
tronic devices 0, El] • On the other hand, in implemen- 
tations such as threshold switches [J, |5j, phase change 
memory Q, and resistive memory Q, CF facilitate in- 
formation storage and logic operations. 

Despite a long history of observations, a theoretical 
framework for CF does not exist. Here, we introduce a 
general thermodynamic theory that describes CF coupled 
with an external circuit and predicts the CF radius as a 
function of the electric current and material parameters, 
as well as the corresponding current-voltage (IV) char- 
acteristics. A finite element numerical simulation is em- 
ployed to support our analytical results, which are found 
to be in excellent agreement with the data. 

For specificity we consider the archetypal CF system of 
chalcogenide glass threshold switches wherein reversible 
switching takes place between highly resistive (amor- 
phous) and conductive (CF) phases. They have recently 
regained interest in connection with 3D stackable phase 
change memory Q. The empirical fact that the CF ra- 
dius increases with current as r a J 1 / 2 has been known 
[3, S] since the seminal work by Ovshinsky [l| , yet it still 
lacks theoretical description. An early approach based 
on the principle of least entropy production |g| did not 
result in verifiable predictions; its validity remains ques- 
tionable [13] and avoidance of it leads to different results, 
as shown here. 

Whether electronic 11, IH or crystalline [HI, or oth- 
erwise different from the host material, CF represents a 
domain of different phase, thus calling upon the analy- 
sis of phase transformations. Our conservative approach 
avoids the principle of least entropy production start- 
ing instead with the kinetic Fokker-Planck equation, the 



applicability of which to phase transformations, partic- 
ularly to nucleation phenomena (according to Zeldovich 
theory) , is well established [14j . Following that approach 
(see e.g. p. 428 in Ref. [I3|) the Fokker-Planck equation 
in the space of cylinder radii r takes the form, 
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Here, / is the distribution function so that f(r)dr gives 
the concentration of filaments in the interval (r, r + dr); 
s is the flux in radii space (s" 1 cm -3 ). The 'fila- 
ment radius diffusion coefficient' B can be estimated as 
va 2 exp(—W a /kT) where v is the characteristic atomic 
frequency (~ 10 13 s" 1 ), a is the characteristic interatomic 
distance, W a is the kinetic phase transformation barrier, 
k is Boltzmann's constant, and T is the temperature. 
A is connected with B by a relationship which follows 
from the fact that s = for the equilibrium distribution 
fo(r) oc exp[— F(r)/kT], where F is the free energy. 

We note that the concept of free energy F that appears 
with the equilibrium distribution /o is not not compro- 
mised by the fact that electric current flows through the 
CF, since that current is fixed by the external circuit and 
serves only as a temperature source. Hence, F describes 
the free energy of the CF in an insulating host paramet- 
rically dependent on the electric current. 

The boundary condition f(r = 0) = to Eq. ([1]) 
reflects the fact that very thin filaments cannot exist due 
to limitations such as loss of conductivity or mechanical 
instability (extraneous to the present model). Another 
condition, f(r = oo) = 0, implies that only finite radii 
are achievable over finite times t. 

Using the right-hand-side expression for s, multiplying 
Eq. ([1} by r, integrating from to oo by parts, and noting 
that J frdr — (r), yields d(r)/dt = {dF/dr). We then 
approximate (F) = F{(r)) and (dF/dr) = d(F)/d(r), 
thereby neglecting fluctuations in the ensemble of nomi- 
nally identical filaments. Omitting for brevity the angu- 
lar brackets, one finally obtains, 
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This equation which expresses the average evolution of 
CF cylinder radius has the standard meaning of a relation 
between the (growth) velocity and the (thermodynamic) 
force —dF/dr, with the mobility b and the diffusion co- 
efficient B obeying the Einstein relation. 

We limit the present discussion to the steady state case 
dr/dt = 0, which, according to Eq. ([2]) takes place when 
the F is a minimum (obviously different from the condi- 
tion of least entropy production 9]). However, in princi- 
ple, Eq. ((2} is capable of describing various transients. 

To analytically present the free energy we consider a 
model in Fig. Q] based on a flat plate capacitor of area A 
and thickness h containing a cylindrical CF of radius r. 
Also, we will assume the characteristic filament dimen- 
sions well above the screening length, thereby neglect- 
ing possible effects due to electric charge redistribution 
around the CF. Including these effects would result in 
additional terms to the free energy adding mathematical 
complexity while not changing the approach of this work. 
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FIG. 1: (color online) Analytical model components with: (a) 
circuit schematic showing source voltage V, load resistance 
Rl, capacitance C, and filament resistance R; and (b) flat- 
plat capacitor of height h, and filament of radius r carrying 
the current /. 

With the above in mind, the major part of the free 
energy is given by, 
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Here, c is the specific heat of the CF, ST is the current de- 
pendent temperature change, U is the voltage across the 
capacitor, C — eA/Anh , e is the dielectric permittivity, 
a is the surface energy, and p is the change in chemical 
potential between the two phases. The first term in Eq. 
dU represents the thermal contribution, the second is the 
electrostatic energy, and the last two correspond to the 
phase transformation. 

For the remainder of this analysis we specify the free 
energy to the case of CF radii r > h because: (I) pub- 
lished experimental results are available; (2) CF dimen- 
sions correspond to the above assumed largeness com- 
pared to the screening length; and (3) the typical values 
of the parameters enable one to neglect the second and 
the third terms in Eq. ([5]), which allows analytical solu- 
tions. In that case, the parameters in Eq. (|3]) are, 
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where V is the source voltage, Rl is the load resistance, 
/ is the current, R and p are the filament resistance and 
resistivity, respectively, and k is the thermal diffusivity, 
taken to be the same for the filament and host materials. 
Substituting Eqs. (@|, the free energy becomes 
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where x = r/ro and we have introduced the dimensionless 
parameters, 

a _ *r 3 V 2 = r_oCV^ R L nr 2 

\2Wnp' h 3W ' ph ' 1 ' 

Here the characteristic energy and length, 

W = 167rcr 3 /3Ai 2 and r = 2a /p. (7) 

would have the physical meaning of nucleation barrier 
and radius in classical nucleation theory (in which p is 
negative and \p\ is used instead). Assuming a and \p\ 
to be of the same order of magnitude as for crystal nu- 
cleation in chalcogenide glasses, one can use the corre- 
sponding estimates [15( W ~ 2 eV and ro ~ 3 nm. 

The free energy as a function of filament radius for 
various source voltages is illustrated in Fig. [H[a). 
The curves indicate that the filament can exist in a 
long-lived metastable state at a; = xt [i.e. with the 
right minimum shallower than the left one at F(x = 
0)]. It becomes stable at relatively high voltage, V > 
(h 2 /CR L )V3Wp/Trr^. On the other hand, finite radius 
filaments become unstable at source voltages below, 



V = m^Wnp/nrl 



(8) 



Vo is defined by the conditions dF/dr — d 2 F/dr 2 = 
and is presented by the curve labeled 0.3 V in Fig. HJa). 
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FIG. 2: (color online) (a) Typical free energy vs. filament 
radius x — r/ro at various source voltages V [from Eq. (5J]. 
The steady state radius is xj and CF are unstable at V < Vo. 
(b) Free energy of a filament as created (thermal contribution 
neglected) and in the steady state. Arrows show how the 
energy minimum moves to the right and becomes metastable, 
separated by a barrier from the state without the filament. 

At V — Vo, the steady state CF radius takes on its 
minimum value, 
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The related filament resistance is a maximum, R max = 
ph/irr^ riin — Rl/2. The characteristic holding current, 
Ih, (below which the CF disappears) can be defined as 
the current attained at the minimum source voltage Vq 
[see the circuit in Fig. Ha)], 
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Along the same lines, the holding voltage is given by 
V h = I h R m ax, which leads to, V h = V /3 (V h is the 
voltage across the bulk of the filament and should not be 
confused with the source voltage Vo). 

The metastable filament is predicted to be extremely 
long-lived at source voltages just slightly above Vo- In- 
deed, as seen from Fig. GJa), the activation barrier 
separating the metastable minimum can be as high as 
Wb ~ 10 8 eV. More quantitatively, close to Vo the shape 
of the free energy is described by the expansion, 



6F = k^— 6xSV 



1 d 3 F 
3! a^" 



\v ,x mi ,S 5x f-, 



where x min = r mm /r , Sx = x - x min and SV = V - V Q 
yielding, 
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The large barrier values are due to a large number of par- 
ticles constituting the filament: h/Hr ~ h> r min/ r 1- 
As illustrated in Fig. &b), the metastable nature of a 
steady state CF does not appear with the filament imme- 
diately upon creation [i.e. when the thermal contribution 
has not yet taken effect and the first term is excluded 
from Eqs. © and ©]. At 
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the stability of the newly created filament is maintained 
solely by the field. While that interpretation implies a 
threshold voltage, it is quantitatively limited by the fact 
that the present model considers filament creation and 
disappearance in one step processes, neglecting the pos- 
sibility of filament nucleation [ItJ • 

Beyond the critical region of voltage close to Vo, the 
filament radius becomes proportional to v7, 
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consistent with the above mentioned experimental obser- 
vations [JIH, (see Fig. EJ). 

Eq. (fT2|) predicts the often observed nearly vertical 
current voltage characteristic where the device voltage, 
Vd = IR oc I/r 2 , remains constant. The value of this 
constant voltage at high current, I ^> i^, is given by, 

V hoo = V /3 3/2 = V h /V3. 
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FIG. 3: (color online) Current as a function of CF radius in- 
dicating the typical I oc r 2 dependence. Excellent agreement 
is obtained between the analytical result of Eq. (|12p . our 
simulation, and data from Ref. @. 
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FIG. 4: (color online) (a) Simulated IV of a filament for the 
device described in Ref. HI at two different load resistances 
Rl- The holding current Ih and voltages Vhoo and Vh are 
shown, (b) Filament current density for pores (active device 
regions) of diameters 25 /im and 50 pm. Saturation occurs 
when the CF fills the pore. The results are in good agreement 
with the data [j| (not shown here). 



The inequality Vhoo < Vh implies a 'knee' in IV char- 
acteristic and a regime of negative differential resistance 
(NDR), as shown in Fig. 0|a). 

The reason for the apparently infinite dynamic con- 
ductivity dl/dV — > oo at V — > Vhoo is that, according 
to Eq. (|I2[) . the CF automatically adjusts its radius to 
maintain a constant current density, 
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It is a bulk material property and independent of current. 
This fact has been known empirically for more than four 
decades. When the CF grows to the device size, satura- 
tion is achieved and the current density increases linearly 
with current [see Fig. E^b)]. 

For the typical parameter values p ~ 0.1 f2-cm, k ~ 
10" 3 cm 2 s-\ h ~ 3000 nm, R L ~ 100 17, e - 10, and 
A ~ 10 10 nm 2 0, Q, B EH , the numerical estimates for the 
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above derived CF holding current Ih ~ 1 mA and cur- 
rent density J ~ 10 A/cm 2 are in excellent agreement 
with the data without any adjusting parameters. On the 
other hand, the predicted holding voltage Vh ~ 0.3 V 
is considerably lower than the measured Vj, ~ 1 V. The 
latter discrepancy could be expected, since our model 
here does not consider blocking electrodes jl, 12| known 
to add a current-independent contribution to the voltage 
across device [H, 0] ■ 

To complement our analytical work, finite element nu- 
merical simulations were performed using the COMSOL 
multiphysics package. The electric field and temperature 
distributions were determined by simultaneously solving 
the coupled current continuity and heat equations, re- 
spectively. The results were then integrated to calculate 
the free energy of Eq. A search algorithm was used 
to determine the minimum free energy in the parameter 
space of applied voltage and CF radius. The simulations 
relied on neither the flat-plate capacitor geometry nor 
the field and temperature approximations of Eq. Q . 

Numerous simulations were performed for a broad 
range of device sizes and geometries. Figs. [3] and U pro- 
vide samples of simulation results for the device structure 
described in Fig. 1 of Ref. d; overall, very good agree- 
ment was obtained without adjusting parameters. The 
simulation results in Fig. EJa) clearly indicate the NDR 
knee in the bottom part of the IV curve, where the CF 
radius increases with I much faster than \fl ' . That re- 
gion of NDR corresponds to the critical region described 
in our analytical treatment [cf. Eq. ([10]>] . 

In summary, we have developed a thermodynamic the- 
ory of steady state CF starting from the basic kinetic ap- 
proach (Fokker-Planck equation) . Analytical expressions 
have been derived and numerical model developed for the 
average CF characteristics. Our results correctly predict 
the filament properties observed in threshold switches 
and the corresponding features of their IV characteris- 
tics, particularly, negative differential resistance vanish- 
ing at high currents; the agreement with experimental 
data is remarkable. Future work will extend this theory 
to include transient analysis and different applications 
(memory devices, shunting in thin-film photovoltaics, di- 
electric breakdown, etc.). 
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